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ABSTRACT 

Aims. We present new medium-resolution high-S/N optical spectra of the ultracompact low-mass X-ray binaries 4U 0614+091 and 4U 1626-67, 
taken with the ESO Very Large Telescope. They are pure emission line spectra and the lines are identified as due to C n-iv and O n-m. 
Methods. Line identification is corroborated by first results from modeling the disk spectra with detailed non-LTE radiation transfer 
calculations. Hydrogen and helium lines are lacking in the observed spectra. 

Results. Our models confirm the deficiency of H and He in the disks. The lack of neon lines suggests an Ne abundance of less than about 10 
percent (by mass), however, this result is uncertain due to possible shortcomings in the model atom. These findings suggest that the donor stars 
are eroded cores of C/O white dwarfs with no excessive neon overabundance. This would contradict earlier claims of Ne enrichment concluded 
from X-ray observations of circumbinary material, which was explained by crystallization and fractionation of the white dwarf core. 
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1. Introduction 

Low-mass X-ray binaries (LMXBs) consist of a neutron star 
or black-hole accretor and a low-mass donor star (M S 1 M Q ). 
Of particular interest are those systems with orbital periods 
Porb 5 80min, which is the minimum period for LMXBs with 
hydrogen-rich main sequence donors. In these ultracompact 
binaries (UCBs) the mass donor must be a non-degenerate 
hydrogen-deficient star or a white dwarf (e.g. Verbunt & van 
den Heuvel 1995). Currently eight such systems with mea- 
sured orbital periods (P rb=l l-50min) are known (Ritter & 
Kolb2003). 

In the recent past, the existence of a group of five ultra- 
compact systems with neon-rich white dwarf donors has been 
claimed based on X-ray spectral properties (Schulz et al. 2001, 
Juett etal. 2001, Juett & Chakrabarty 2003). The donors are 
then C/O or even O/Ne/Mg white dwarfs that have transferred 
a significant fraction of their mass to the neutron star, in con- 
trast to the usual wisdom that the donors are the remains of He 
white dwarfs. This has caused new explorations of the forma- 
tion of these systems (e.g. Yungelson etal. 2002). Our motiva- 
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tion for studying these systems is that the stripped donor stars 
offer the possibility to probe the interior composition of white 
dwarfs, which depends on the interplay of gravitational settling 
and crystallization of chemical elements. 

Three of these five Ne-rich systems belong to the above- 
mentioned eight UCBs with measured orbital periods, while 
two of them are believed to be UCBs because of their similar 
optical and X-ray properties. The class of Ne-rich UCBs con- 
sists of: 

(1) 4U 1626-67 (P orb =41min) 

(2) 4U 06 14+091 

(3) 2S 0918-549 

(4) 4U 1543-624 (P orb =18min, Wang & Chakrabarty 2004) 

(5) 4U 1850-087 (P orb =20min) 

And 4U 1626-67 may be regarded as the prototype of this class. 
The donor's Ne-rich C/O- WD nature is derived from X-ray and 
ultraviolet spectra that exhibit double-peaked emission lines 
that obviously stem from the accretion disk (Schulz et al. 2001, 
Homer etal. 2002). 

The close relation of the other four objects (2) -(5) to 
4U 1626-67 was based on extraordinary high Ne/O abundance 
ratios (when compared to the ISM value) measured from 
ASCA spectra of (3)-(5) and a Chandra spectrum of (2). The 
spectra exhibit photoelectric absorption edges of neutral O and 
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Fig. 1. Blue spectra of 4U 0614+091 and 4U 1626-67. H and He lines are lacking. The emission lines are identified as Oii-iii 
and C ii-iii (see also Table|2}. All absorption features are of interstellar (i.s.) origin: Ca n 3934/3968 A, the 4300 A CH band, and 
two diffuse interstellar bands (DIB). Plotted over the two observed spectra is a synthetic accretion disk spectrum reddened with 
E(B-V)=0.55 and E(B-V)=0.40, respectively. 



Ne in the interstellar medium (ISM) along the line-of-sight, 
which is suspected to originate from expelled material close 
to the binary systems. New X-ray spectroscopic observations 
of (3) and (4) obtained with Chandra and XMM-Newton con- 
firmed the ASCA results, although different values for the 
Ne/O enrichment have been obtained from a Chandra and a 
XMM-Newton spectrum of (4) (Juett & Chakrabarty 2003). 
In contrast to the ASCA measurement of (5), recent XMM- 
Newton and Chandra observations found no evidence of an un- 
usual Ne/O ratio (Sidoli et al. 2005, Juett & Chakrabarty 2005). 
These apparently contradictory results can be attributed to a 
variable Ne i/O i ratio due to changes in the ionisation structure 
in the measured absorption columns that, however, are not un- 
derstood (Juett & Chakrabarty 2005). Hence, this means that 
the measured Ne/O ratio does not reflect the donor composi- 
tion. Although the orbital periods of these systems are below 
80 min (or at least believed to be so small), this does not nec- 
essarily mean that they contain C/O donors. For example, the 
X-ray burst properties of the 4U 1820-30 (P or b=l 1 min) suggest 
an He-WD donor in that ultracompact system (e.g. Strohmayer 
& Brown 2002). Another example are the relatively short X-ray 
bursts observed from (3) that even suggest H-burning of ma- 
terial accreted onto the neutron star (NS) (Jonker etal. 2001), 
while optical spectra suggest a C/O WD donor (Nelemans et al. 
2004 and this work). A possible explanation is that the heavy 
elements (C, O, Ne) undergo spallation during accretion leav- 
ing H and He nuclei for thermonuclear burning on the NS 
(Bildsten etal. 1992). In a recent paper In't Zand etal. (2005) 
conclude that the Ne/O ratio and the X-ray burst properties are 
all best explained with an He-rich donor. 



In conclusion, the only way to confirm that the four sys- 
tems (2)-(5) - besides the prototype (1) - indeed contain C/O- 
rich donors, perhaps enriched with Ne, is by UV and/or opti- 
cal spectroscopy. Nelemans et al. (2004) show that the optical 
spectra of (2)-(4) are devoid of hydrogen and helium emis- 
sion lines, which are usually seen in H- or He-rich accreting 
systems. Their spectra exhibit low-ionisation C and O emis- 
sion lines, most prominent in the brightest of these objects, 
4U 06 14+091. In this paper we present new optical spectra of 
this system that cover a larger wavelength interval. They con- 
firm the earlier conclusion by Nelemans et al. (2004) that the 
emission lines arise from the C/O-dominated material that is 
probably located in the accretion disk and not in the X-ray 
heated face of the white-dwarf donor. We also present a first 
detailed optical spectrum of the prototype 4U 1626-67, which 
mainly shows weak oxygen emission lines and proves the H- 
and He- deficiency in this system, too. Nelemans & Jonker 
(2005) also performed VLT observations of this system (also 
in spring 2004) with a similar setup. They show a section of 
their spectrum and emphasize the similarity with 4U 06 14+91. 

In addition we present results from first attempts to model 
the observed spectra with non-LTE accretion-disk models. We 
derive upper limits for the H and He abundances and inves- 
tigate the formation of neon lines. Detailed C and O line- 
formation calculations can already qualitatively explain the ob- 
served emission lines. The ultimate goal is to derive detailed 
abundances and other disk parameters from the observed line 
profiles. 

The paper is organised as follows. We describe our observa- 
tions in the following section. We then present our line identi- 
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Fig. 2. Red spectra; 4U 0614+091 shows emission lines from Cii-iv and Oii-m, while the 4U 1626-67 spectrum is virtually 
featureless. All absorption features are either telluric (marked by horizontal bars) or of interstellar origin (Nai 5890/5896 A and 
a DIB at 5780 A). The features at 5575 A are artifacts ("x") due to an O i sky emission line. Overplotted is the synthetic accretion 
disk spectrum like in Fig.^ 



fications in Sect. 3. Section 4 contains a description of our disk 
model calculations and we summarise and conclude in Sect. 5. 



2. Observations 

We obtained optical medium-resolution long-slit spectra of 
4U 1626-67 and 4U 0614+091. The V magnitude of both sys- 
tems is 18.5 (Ritter & Kolb 2003). We used the FORS1 spectro- 
graph attached to UT1 of ESO's Very Large Telescope (VLT) 
on Paranal in Chile. The slit width was 1". We used two grisms 
(600B and 600R, the latter in combination with order separa- 
tion filter GG435) and obtained spectra covering the regions 
3600-6000 A and 5400-7500 A with a mean dispersion of 1 .20 
and 1.07 A pix~', respectively. Observations were performed 
in service mode between Nov. 2003 and Mar. 2004. Each tar- 
get was observed at least twice, see Table [l\ for details. For 
4U 1626-67 the exposure time covers almost one orbital pe- 
riod so that any respective spectral variation is smeared out. 
The data were processed through the standard ESO reduction 
pipeline. In addition we performed a flux calibration using ob- 
servations of the DA white dwarf EG 274 with the same in- 
strumental setup. Since this flux standard was observed only 
once per wavelength range and since the science targets were 
mostly observed at different dates, this provides only a rough 
absolute flux calibration. For each object the spectra were co- 
added to obtain one final spectrum. The blue and red spectra 
are displayed in Figs.^andlJ] respectively. 

Archival spectra of 4U 1626-67 taken with the Hubble 
Space Telescope (HST) and the STIS instrument cover the 



complete UV/optical wavelength range (1150-10000 A). The 
FUV spectrum is described in detail by Homer etal. (2002) 
and a first comparison with synthetic accretion disk spectra 
was presented by Werner etal. (2004). In Fig.[3]we compare 
the optical HST spectrum taken from the MAST archive with 
our VLT spectrum on an absolute flux scale. The blue and red 
sections of the VLT spectrum were scaled by a factor of 1 .65 
and 1.45, respectively, to normalise them to the HST flux at 
5600 A. It is obvious that the VLT spectrum is not only weaker 
but also flatter than the HST spectrum. This can be attributed 
either to flux calibration problems or to source variability. 

3. Line identification 

Emission lines from hydrogen and helium are completely lack- 
ing in the VLT spectra of both systems. We detect neither 
Balmer lines nor lines from He i (e.g. 4471 or 5876 A) or He n 
(e.g. 4686 A), which are usually seen in optical spectra of in- 
teracting binaries. In the following two sections we discuss the 
line spectra of both systems in detail. Line identification is per- 
formed "by eye" without an equivalent-width estimation. 

3.1. 4U 0614+091 

The spectrum of 4U 0614+091 shows numerous emission fea- 
tures that can be assigned to ionized carbon and oxygen, 
namely C ii-iv and O n-m. The detected emission features are 
listed in Tabled together with possible line identifications. 
Most features are blends of lines from at least two ions, making 



4 K. Werner et al.: VLT spectroscopy of LMXBs 4U 0614+091 and 4U 1626-67 

Table 2. Emission features observed in the spectra of 4U 1626-67 and 4U 0614+091 and suggested line identifications. A plus 
sign in parentheses denotes an uncertain detection. 
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Fig. 3. Complete VLT spectrum of 4U 1626-67 compared to a spectrum taken with HST. The VLT spectrum, scaled to the HST 
spectrum at 5600 A, is flatter possibly due to problems in absolute flux calibration. The HST spectrum has a poorer resolution 
and S/N-ratio, although some emission features at A < 5000 A can be recognized in both datasets. Artificial emission spikes are 
marked by "x". Overplotted is the synthetic accretion disk spectrum, like in previous figures, but now reddened with E(B-V)=0.20 
to fit the continuum shape of the HST spectrum. 
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Fig. 4. Detail from three single spectra of 4U 0614+091. Thick line: Spectrum taken on Nov. 03, 2003. Two thin lines: Spectra 
taken consecutively on Dec. 15, 2003 (see Table The strength of the two strongest emission features, the Cn multiplets at 
6580 A and 7235 A, show an increase over a time interval of several weeks. 



interpretation quite difficult, but some of them are probably un- 
blended. These are the features at 5590 A (O m), 5650 A (C n), 
5695 A (C m), and 7235 A (C n). The emission features at 4710 
and 4940 A are perhaps exclusively due to O n, but they are 
blends of four and two multiplets from this ion, respectively. 
Nelemans et al. (2004) present optical spectra of this star taken 
with VLT+FORS2 with slightly different wavelength coverage 
(4620-8620 A). We essentially confirm the detection of their 
emission features and line identifications, however, our spectra 
extend to shorter wavelengths, down to 3600 A, which allows 
us to detect some additional emission lines. On the other hand, 
we do not see four of the weak emission features identified in 
the Nelemans et al. (2004) spectra (at AA 5 140, 5 190, 5280, and 
6070 A). This could be the consequence of line variability. In 
Fig.0]we compare our single spectra of 4U 0614+091 in detail, 



and it is possible that the emission strength of the strongest C n 
lines is variable on a time scale of weeks. In order to assess 
the significance of this variability, we simulated several times 
with different noise the strongest emission features observed 
on Nov. 3 and Dec. 15, applying the S/N ratio of the VLT spec- 
tra. We find a 2<x probability that the variability is real. The 
line widths correspond to a projected rotational velocity of the 
order 1200 kms~'. 

3.2. 4U 1626-67 

The spectrum of 4U 1626-67 clearly shows less features than 
does that of 4U 06 14+091. Only five emission lines can be 
detected in the blue spectrum (see Tabled} and they are seen 
in both binaries. The red spectrum of 4U 1626-67 is virtually 
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Table 1. Log of the observations. 



Date, UT 


Grism 


T 

1 exp 


Airmass 


Seeing 






(sec) 




(arcsec) 




4U 0614+091 






2UUJ Nov. U4, U/:2z:Uo 


600B 


1735 


1.22 


0.73 


ZUUJ 1\UV. '-'^5 U / . Jt 1 . J / 


600B 


1735 


1.20 


78 
u. / o 


2003 Nov. 03, 07:40: 17 


600R 


1735 


1.21 


1.00 


2003 Dec. 15, 06:25:42 


600R 


1735 


1.28 


0.52 


2003 Dec. 15, 07:00:31 


600R 


1735 


1.37 


0.57 




4U 1626-67 






2004 Mar. 23, 08:22:00 


600B 


1735 


1.37 


0.76 


2004 Mar. 23, 09:01:42 


600B 


1735 


1.36 


0.77 


2004 Mar. 21,08:48:53 


600R 


1735 


1.36 


1.58 


2004 Mar. 22, 08:19:40 


600R 


1735 


1.37 


0.87 


2004 Mar. 22, 09:00:00 


600R 


1735 


1.36 


0.71 


2004 Mar. 23, 07:41:51 


600R 


1735 


1.40 


0.57 




EG 274 






2004 Mar. 23, 09:48:17 


600B 


50 


1.05 


0.67 


2004 Mar. 23, 09:53:25 


600R 


40 


1.05 


0.67 



continuous. The strong Cn emissions seen in 4U 0614+091 at 
6580 A and 7235 A are lacking in 4U 1626-67. That could point 
to a higher ionisation of carbon in 4U 1626-67, although the 
Civ 5810 A line is lacking, too. The presence of On lines in 
the blue spectrum of 4U 1626-67 (e.g. at 4345 A) also contra- 
dicts a higher ionisation, rather, it appears that the differences 
in the spectra of both binaries could be assigned to a differ- 
ent C/O ratio, being higher in 4U 06 14+091. The entire set of 
emission lines in 4U 1626-67 could be due to oxygen alone, but 
in the UV spectrum C is clearly present (Homer et al. 2002). In 
addition, the line widths appear broader than in 4U 0614+091, 
however, this is difficult to quantify because of possible line 
blends and an uncertain identification of the continuum. 

All absorption features seen in our optical spectra are ei- 
ther of interstellar origin or telluric. They are marked in Figs. [2 
and|2 

4. Exploratory disk models 

We began the construction of accretion disk models to calcu- 
late synthetic spectra and report here on the current state of our 
work. We use our accretion disk code AcDc, which is described 
in detail by Nagel etal. (2004). In essence, it assumes a radial 
c-disk structure (Shakura & Sunyaev 1973). Then the disk is 
divided into concentric annuli. For each annulus we solve the 
radiation transfer equation (assuming plane-parallel geometry) 
together with the non-LTE rate equations for the atomic level 
populations, plus energy- and hydrostatic equations, in order to 
calculate a detailed vertical structure. The integrated disk spec- 
trum is then obtained by co-adding the specific intensities from 
the individual annuli, accounting for inclination and Keplerian 
rotation. 

It is not the aim of this paper to present a detailed fit to the 
emission line spectra of the two binaries. This requires exten- 
sive parameter studies that are extremely time consuming. We 
rather choose to calculate synthetic spectra from selected disk 



Table 3. Summary of model atoms used in the disk model cal- 
culations. For each ion we list the number of NLTE levels and 
the number of line transitions. In brackets we give the number 
of lines in the wavelength range covered by our optical spec- 
tra after fine-structure splitting for detailed line-profile calcula- 
tions. 



Element 


Ion 


NLTE levels 


Lines 


H 


1 


10 


45 (8) 




ll 


1 
1 




He 


I 


29 


61 (16) 




ll 


14 


91 (14) 




ill 


1 




C 


I 


7 


4(0) 




ll 


38 


160 (42) 




in 


58 


329 (143) 




IV 


9 


17(2) 




V 


1 




O 


I 


1 






11 


29 


82 (56) 




III 


36 


42 (47) 




IV 


11 


5(0) 




V 


6 


4(0) 




VI 


1 




Nc 


I 


3 







11 


68 


232 (49) 




III 


4 


0(0) 




IV 


1 





regions in which we believe the physical properties are repre- 
sentative for the formation of the observed optical spectra. This 
should at least give a rough idea of the relative strength of lines 
from different ionisation stages and, thus, is primarily thought 
to put our line identifications on firm ground. At present we 
neglect irradiation of the disk by the neutron star, because it 
would introduce new free parameters. We expect that the ra- 
dial ionisation structure of the disk will be shifted to larger 
radii when irradiation is taken into account, but we hope that 
the relative line strengths are not affected to the extent that our 
identifications become completely wrong. Although the stud- 
ied systems are strong X-ray sources, we do not see recombi- 
nation lines from highly ionised species in the optical spectra 
as might be expected. It is conceivable that the outer parts of 
the disks in which the optical spectrum arises is shielded from 
X-ray irradiation by an inflated inner-disk region. 

Since we do not compute the spectrum of the entire disk, 
we cannot expect to match the overall observed continuum flux, 
and we have already pointed out that there are problems with 
the absolute flux calibration. We use interstellar reddening as a 
free parameter in order to roughly fit the model to the observed 
flux level. The applied reddening is given in the figure captions. 

The synthetic spectra presented here are based on the fol- 
lowing choice of disk-model parameters. The central object is a 
neutron star with 10 km radius and a mass of 1 .4 M . The mass 
accretion rate is 2 • 10~ 10 M Q /yr. The disk extends from 1000 
to 2000 stellar radii. The corresponding Keplerian velocities at 
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Fig. 5. Vertical structure of the C/O-disk model at a distance of 
20 000 km from the neutron star. The emergent disk flux at this 
location corresponds to T e ff =28 000 K. The physical variables 
are plotted against the column mass measured from the surface 
towards the mid-plane. 



these radii amount to 4300 and 3000 kms" 1 , and the effective 
temperatures to 47 000 K and 28 000 K, respectively. The disk 
is divided into five annuli such that T e fi decreases almost lin- 
early with radius steps. The inclination angle is set to 10°. The 
Reynolds number used to parametrise the disk viscosity was set 
to 10000, which corresponds to a sa 0.5. The chemical compo- 
sition is C=50% and O=50% (by mass). 




-5 -4 -3 -2 -1 

log (m / g cm' 2 ) 



Fig. 6. Vertical ionisation stratification of chemical elements in 
the disk models at a distance of 20000 km from the neutron 
star. The vertical line drawn at log m = -0.9 indicates tr oss = 1 . 

For the opacity and emissivity calculations, it is essential to 
solve the non-LTE rate equations with detailed model atoms. 
Our main emphasis was put on the C ii-m and O ii-m ions. The 
number of non-LTE levels and radiative line transitions are 
summarised in Table|3] Level energies, oscillator strengths, and 
bound-free cross-sections for photoionisation are taken from 
the Opacity Project (Seaton etal. 1994) TOPbase 1 . Electron 
collisional rates for (de-) excitation, ionisation, and recom- 
bination are computed with the usual approximate formulae. 

1 http://vizier.u-strasbg.fr/topbase/topbase.html 
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Fig. 7. Comparison of models with zero and 10% hydrogen and 
helium content to the observations. From the lack of H^g and 
He n emission lines in the observed spectra, we conclude that 
the accretion disks are strongly H and He deficient. The model 
spectra were normalised to the local continuum flux. 
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Fig. 8. Comparison of models with zero and 10% neon content 
to the observations. The model including neon exhibits weak 
Ne ii lines in the region 4350 - 4450 A, which cannot be de- 
tected in the observed spectra. 



The neon model atom was taken from Dreizler (1993). For 
the final spectrum synthesis, fine structure splitting of atomic 
levels must be considered; level populations were taken from 
the models with appropriate statistical weights. Level energies 
were obtained from the NIST 2 database. The optical synthetic 
spectra finally include a total number of 377 lines from Hi, 

2 http://physics.nist.gov/ 



He i-n, C n-iv, O ii-m, and Ne n. They essentially contain most 
of the possible C and O line identifications given in Table |5J 
plus many others that turn out to be too weak to be seen in the 
rotationally broadened spectra. 

As an example we present the vertical structure of the ac- 
cretion disk model at a distance of 2000 stellar radii from the 
neutron star. The emergent flux at this location corresponds 
to T e s =28 000 K. Figure [5] shows the run of several quanti- 
ties above the disk midplane on a column mass scale m. The 
Rosseland optical depth reaches unity at login = -0.9. Here 
the gravity amounts to logg=5.8. Mass and electron densities 
are 6 • 10 8 gem" 3 and 5 • 10 15 cm 3 , respectively. Hence the 
physical conditions in the line-forming regions are comparable 
to those in the atmosphere of a subluminous B star. The verti- 
cal distribution of ionic fractions of C and O is shown in Fig. [6] 
There we also show the H, He, and Ne fractions that were taken 
from the test models that include these species in an amount of 
10%. The dominant ionisation stages of oxygen at tr oss = 1 
are On and Om, giving rise to prominent emission lines of 
these ions in the spectra. In the case of carbon, C m dominates, 
closely followed by C n and C iv. This explains why we see 
lines from three C ions. The dominant neon ions in the line- 
forming regions are Ne n and Ne in. For helium we have He n 
dominating, followed by He i, so we would expect strong lines 
from these ions if neon and helium were abundant. Hydrogen 
is ionised by about 99.9%, but still, prominent emission lines 
are predicted if H were an abundant species (see below). 

4.1. C and O lines compared to observations 

The synthetic spectrum is plotted together with the observed 
spectra in Figs. [2 and [2 Let us first consider 4U 06 14+091. 
Generally, many of the observed features are also seen in the 
model, although they do not match the strength. This basically 
corroborates our line identifications. The C n lines of the model 
are too weak (e.g. at 6580 A and 7235 A), while the C in emis- 
sions are in some cases too weak (e.g. 4650 A) or in other cases 
too strong (4180 A). The Civ 5810 A line matches well. The 
oxygen lines show a similar behaviour. The On lines of the 
model are too weak (e.g. at 4940 A). Some lines of O m match 
reasonably well (e.g. at 3720 A), while others do not (e.g. at 
5590 A, which is much too weak in the model). 

The model comparison with 4U 1626-67 confirms our ideas 
from the first inspection of the spectrum. The lines in the 
blue region can be explained qualitatively by the mere pres- 
ence of oxygen lines. This, and the missing Civ 5810 A line 
might originate in a relatively low C/O ratio when compared to 
4U 06 14+091. 

It is disappointing that our NLTE disk model obviously 
gives a poorer fit to the observed line spectrum than the sim- 
ple isothermal, constant-density, LTE slab model presented by 
Nelemans et al. (2004). However, that slab model can at best 
fit a a limited spectral range, as it is emitted from a partic- 
ular emission region with an assumed value for temperature 
and density and will never be able to simultaneously fit the ob- 
served spectra from the X-ray to the UV and optical ranges. In 
good agreement with the observation, our model exhibits lines 
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from three ionisation stages of carbon, indicating that temper- 
ature and density in the line-forming regions are reproduced 
reasonably well. The fact that the strength of many emission 
lines of a particular ion is either over- or underestimated may 
stem from drawbacks in the model atoms. One reason could be 
that the model ions are still too small and ignore the interlock- 
ing effects of neglected energy levels. A further extension of 
the model ions is hampered by the lack of atomic data, mostly 
oscillator strengths. Another reason for the poor line fits could 
be errors in the electron collisional rates. Only few are known 
from experiments or quantum mechanical calculations. 

These problems also might affect our neon line-formation 
calculations. However, the employed Ne model atom has been 
designed for the analysis of sdO stars and gave successful fits 
to observed Ne n lines (Dreizler 1993). Therefore we think that 
the predicted Ne emission lines in the model are more reliable 
than the majority of the C and O lines. Atomic data for H and 
He are accurately known and our model atoms have been em- 
ployed with success to analyse many classes of stellar spectra, 
hence, we regard the computed line strengths of H and He as 
much more reliable than those of the metal lines. Another rea- 
son for the poor line fits might be that the assumed a-disk does 
not describe the physics of the emitting region well. 

4.2. Limits on abundances of hydrogen, helium, and 
neon 

The upper abundance limits we derive for H and He are re- 
garded as realistic because of the above considerations, but the 
limit for neon is less secure. 

In addition to the pure C/O disk model we also computed 
a model including 10% hydrogen and a model including 10% 
helium in order to see if this allows the H and He abundances to 
be constrained. In Fig.^we display the resulting emission lines 
of Hg and Hen 5412 A. From the lack of observed emission 
lines, it is evident that H and He are at most present at the 10% 
level. This excludes the possibility that the disks of 4U 1626- 
67 and 4U 06 14+09 1 are dominated by hydrogen or helium and 
confirms the suspicion that they are in fact C/O-dominated. 

We also computed a C/O model that includes 10% neon. 
It exhibits weak Ne n lines, but no such lines are detected in 
the VLT spectra (Fig. [3}. We conclude that the neon abundance 
cannot be larger than « 10%. 

5. Summary 

We have presented new high-quality optical spectra of the ultra- 
compact low-mass X-ray binaries 4U 0614+091 and 4U 1626- 
67. They are pure emission line spectra and most probably stem 
from the accretion disk. The spectral lines are identified as due 
to C n-rv and O n-nx Line identifications are corroborated by 
first results from modeling the disk spectra with detailed non- 
LTE radiation transfer calculations. Hydrogen and helium lines 
are lacking and our models confirm the deficiency of H and He 
in the disk. Hence, the donor stars in these systems are in fact 
the eroded cores of C/O white dwarfs. There are indications 
that the O/C ratio in 4U 0614+091 is higher than in 4U 1626- 



67. This could suggest that the stripping process of the WD in 
4U 06 14+091 is more advanced. 

It is hard to estimate the systematic error of the derived up- 
per limit for the neon abundance so that the following conclu- 
sions are at the moment rather uncertain. From the lack of Ne n 
lines we find that the Ne abundance cannot exceed w 10% or, 
in other words, the Ne/O ratio is at most of the order 0.2. The 
much higher Ne/O ratios (w 0.7) derived from the ISM X-ray 
absorption edges of neutral Ne and O (Schulz etal. 2001, Juett 
etal. 2001) would produce detectable Nen lines in the disk 
spectra. This confirms the conclusion of Juett & Chakrabarty 
(2005), that the determined ISM abundances of Ne and O are 
affected by ionisation effects and, hence, do not reflect the 
abundances of the donor stars. 

For an initial solar metallicity the 22 Ne abundance in an 
Ne-enriched crystallized and fractionated WD core can be esti- 
mated from theoretical models to m 0.07, but this value is very 
uncertain (Yungelson et al. 2002). It could be even higher by a 
factor of 3 (Isern et al. 1991). Such an Ne-rich core would have 
a mass of ss 0.06 M Q (Yungelson et al. 2002). The mass of the 
WD donor in 4U 1626-67 is much smaller (0.01 M Q ; Yungelson 
et al. 2002). If we accept an upper limit of Ne=0.1 for the ob- 
served neon abundance, then we may draw the following con- 
clusion. Either the WD core has crystallized and fractionated, 
then our observation favors a relatively small Ne-enrichment 
as a result of the crystallization process. Or, if one accepts that 
fractionation would result in a high Ne abundance of 0.2 (unob- 
served), then the WD core in 4U 1626-67 had no time to crys- 
tallize which, depending on various details, lasts several Gyr 
(e.g. Hernanz etal. 1994). 

Future work will concentrate on the disk modeling for the 
prototype 4U 1626-67, for which the observational database is 
the best of all such systems. Ultimately, the spectral properties 
(flux distribution and emission line strengths) over the com- 
plete wavelength range comprising spectral observations with 
Chandra, HST, and VLT, must be explained by a unique disk 
model. 
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